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Abstract 

The infinite source Poisson arrival model with heavy-tailed workload distributions 
has attracted much attention, especially in the modeling of data packet traffic in com- 
munication networks. In particular, it is well known that under suitable assumptions 
on the source arrival rate, the centered and scaled cumulative workload process for 
the underlying processing system can be approximated by fractional Brownian mo- 
tion. In many applications one is interested in the stabilization of the work inflow to 
the system by modifying the net input rate, using an appropriate admission control 
policy. In this work we study a natural family of admission control policies which keep 
the associated scaled cumulative workload asymptotically close to a pre-specified lin- 
ear trajectory, uniformly over time. Under such admission control policies and with 
natural assumptions on arrival distributions, suitably scaled and centered cumulative 
workload processes are shown to converge weakly in the path space to the solution 
of a d-dimensional stochastic differential equation (SDE) driven by a Gaussian pro- 
cess. It is shown that the admission control policy achieves moment stabilization in 
that the second moment of the solution to the SDE (averaged over the d-stations) 
is bounded uniformly for all times. In one special case of control policies, as time 
approaches infinity, we obtain a fractional version of a stationary Ornstein-Uhlenbeck 
process that is driven by fractional Brownian motion with Hurst parameter H > i. 
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tribution, self-similarity, fractional Brownian motion, fractional Ornstein-Uhlenbeck 
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1 Introduction 

This work is motivated by a general workload model for data traffic in communication 
networks considered by Kurtz [6j. In this model, a large number of sources input work into 
a system. Let N(t) be the number of source activations up to time t. For the activation 
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of the ith source, let Xi(s) denote the cumulative work input into the system during the 
first s units of time that the source is on, and let T{ denote the length of time that the ith 
source remains active (ith session length). Then, the total work input into the system up 
to time t is given by 

W(t)= f X n{s) (t n(s) A(t-s))dN(s). (1.1) 

J 

It is assumed that pQ,Tj) are i.i.d. with distribution v on .Dr + [0,oo) x [0,oo), where 
Dr + [0, oo) denotes the space of R + - valued right-continuous functions on [0, oo) having 
left limits (RCLL). 

The paper [B] assumes that the arrival (counting) process N(t) is characterized by an 
intensity A, which may depend on the past of N, W and another process Q, representing 
the number of active sources. For simplicity and since only this case will be considered 
below, assume that the intensity A depends only on the total workload W(t) and time t 
as 

X = X(t,W(t)), (1.2) 

and focus only on the process (N, W). Letting £ be a Poisson random measure on [0, oo) x 
D]r[0,oo) x [0,oo) with intensity measure r] = m x v (m being Lebesgue measure), the 
process (./V, W) can be represented as 

N(t) = a(B(t)), 



W(t) = u(r A(t--y(s)))£(ds,du,dr), (1.3) 

JB(t) 

where ^ 

B(t) = {(s,u,r) : s < [ X(z, W(z))dz} (1.4) 



and 7(5) satisfies 



7( X(z,W(z))dz = s. (1.5) 



See Section [2TT1 for some discussion on this system of equations. The key results of [6] are 
the law of large numbers and the central limit theorem for the scaled system (X n ,Y n ), 
where 

X n (t) = -N n (t), Y n (t) = -W n (t), (1.6) 
n n 

and (N n ,W n ) are defined as in (j 1 . 3 1) but using a Poisson random measure £ n with in- 
tensity measure nm x v and an intensity X n (t,w(t)) = nX(t, ra -1 u;(£)). Under suitable 
assumptions, (X n ,Y n ) converges in probability to (X,Y) satisfying 

X(t)= I X(s,Y(s))ds,Y(t) = f fi{t-r)X(s,Y(s))ds, (1.7) 



Jo 



where fi(t) = E(Xj(ri A t)) (Theorem 2.1 in |6j). Under suitable assumptions, the scaled 
and centered process (X n ,Y n ) = ^fn{X n — X, Y n — Y) converges in distribution to (X, Y) 
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satisfying 

X(t) = E(fl(t)) + f X y (s, Y(s))Y(s)ds, 
Jo 

Y(t)= u(r A(t-j{s)))E(ds,du,dr) + y.{t - s)X y {s,Y(s))Y(s)ds, (1.8) 

J B{t) Jo 

where B(t) and 7(s) are defined as in (|1.4p - (jl.5p but replacing W by Y, and H is a 
Gaussian random measure with the control measure 

E(\E(ds,du,dr)\ 2 ) = dsv{du,dr) (1.9) 

(Theorem 2.2 in [6]). 

One special case of (ATj,Tj) is particularly interesting in the context of modeling data 
traffic in modern communication networks. This is the case where 

Xi(s) = s (1.10) 

and Tj are heavy tailed with the distribution 

u (dr) = {13- l)0(0r + l)~^dr, r > 0, (1.11) 

where /3 6 (2, 3) is the tail index and G (0, oo) is a scale parameter. (Note that we use 
the same notation v to denote the distribution of Tj alone.) When intensity A is constant 
and (|l.lU|) - (ll.lip are assumed, it is well known that, under suitable assumptions and 
proper scaling, the cumulative workload process converges to fractional Brownian motion. 
A version of this fact appears in Section 4 of [6j. But see also [3J 01 O [U HO] . 

For the convergence to fractional Brownian motion, in a scaled system, it is also 
necessary to rescale the measure v(dr) in (II. lip . One way to see this is to observe that 
without rescaling, the Gaussian random measure S in (jl.8p - (|1.9p is not self-similar in the 
variable r. For this reason (see for example, [6], Section 4 for the case when A is constant), 
it is natural to scale the measures as: 

v n {dr) = (/3- I)n6(n6r + ly^dr. (1.12) 

This case is not included in Theorems 2.1 and 2.2 of [6] which, although allow for state 
dependent A, treat the scaled system (jl.6p that has no scaling in the intensity measure v, 
and hence the key convergence results of [6] for non-constant A cannot be applied with 
(|1.12p . In fact, as already suggested by the result in Section 4 of [6] for the constant A 
case, dealing with (|1.12p for non-constant A is expected to be more involved. For example, 
a natural normalization in this case is no longer yfn. 

Models where the arrival intensity A is a function of the state process are natural 
when one considers control mechanisms for regulating the amount of work in the system. 
A common form of a control policy that aims to appropriately balance long processing 
delays with low processor utilization, consists of suitably decreasing the input rate when 
the workload in the system is very high and increasing the rate when it drops too low. 
Study of asymptotic behavior of the workload process with heavy-tailed session length 
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distributions, under such state feedback control mechanisms is the subject of current 
work. We shall consider a scaled multidimensional system where the session lengths are 
distributed according to v n as in (|1.12p . and establish limit theorems for settings where 
A is state dependent. We are particularly interested in the design of control policies 
that keep the net workload (asymptotically) close to a pre-specified linear trajectory 
such that the variability (suitably scaled) is bounded uniformly in time. The slope of 
the linear trajectory represents the system processing rate and thus such control policies 
yield uniform in time reliability bounds on probabilities of processor underutilization and 
overload. 

Let us now describe briefly our model and results established below. We suppose 
that a system consists of d processing stations, and that workload arrives to each station 
(independently of others) as before. The intensity A controlling the arrival rate, however, 
now depends on the average total workload across all the stations. More specifically, 
denoting the total cumulative workload at the ith station by Ui(t) and their average y(t), 
we suppose that 

X = f(t,y(t)). (1-13) 
More specifically yet, we will work with special intensities A having the form 

A = f(t, y(t)) = exp{-g(y(t) - bt)}, (1.14) 

for some b > and function g. The constant b represents the processing rate at each 
station, although processing of work is not explicitly included in our model and plays no 
role in the analysis. The function g will satisfy the assumption stated next. 

Assumption 1.1 g(0) = 0. The function g is twice differentiable and its first and second 
derivatives g' and g" satisfy 

0<£< g'(x) <L, Vx G R, 

and 

\g"(x)\ < L, Vx g R 

for some i, L G (0, oo) . 

The above assumption will be taken to hold throughout this work and will not be explicitly 
noted in the statements of various results. Note that under this assumption g is a strictly 
increasing function and g(u) > if u > and g(u) < if u < 0. From the properties of g, 
we see that the intensity A in (|1.14p has a natural physical interpretation: the intensity 
of session arrivals at the ith station increases when y(t) drops below bt while it decreases 
when y(t) exceeds bt. We will refer to g as an admission control policy. 

Our scaled system will be characterized by independent Poisson random measures 
having common intensity measure n a m x u n where v n is as in (j!.12[) and the cumulative 
workload process Y n ^{t\ unlike (jl.6p will now be normalized by a factor of n Q_1 , rather 
than n (see (|2,6p ). We will assume that 

a£ (/3-l,min{3/3-5,5-/3}). (1.15) 
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The reason for such choice of a, and for the normalization n a 1 will be given below (see 
Remark E2]). 

Precise evolution equations for Y n j are given in Section [2j We now give a brief de- 
scription of our main results. In Theorem 12.21 we prove a law of large numbers result 
stating that, as n — y oo, Y n — (^n,i> • • • >Y niC i) converges in probability, in D^d [0,oo) 

to a continuous (non-random) trajectory U = (U, . . . , U) T , where U is characterized as 
the unique solution of an ordinary differential equation (ODE) (see (|2.8j) ). and a rate of 
convergence is given as well. The solution U has the property that sup t>0 \ U(t) — bt\ < oo. 
In fact, with a particular choice of b, namely b = gm^ m > w e nave U(t) = bt for all t. 

Next, we study the fluctuations of Y n . In Theorem l2.3l we show that suitably centered 
and normalized form of Y n , denoted as Z n (see (12. lip ), converges in distribution, in 
L* K d[0,oo) to the solution Z = {Z\, . . . , Zd) T of a <i-dimensional stochastic differential 
equation (see (I2.12P ). driven by d independent Gaussian processes Ri, i = 1, . . . , d. The 
moment stabilization property of the admission control policy is demonstrated in Theorem 
12.41 which says that sup t>0 E|Z(i)| 2 < oo, where Z = 4 Y2i=i 

We remark that, in the case when b = i one can achieve the law of large number 

limit of bt by simply taking the admission control policy to be g = (this function obvi- 
ously does not satisfy Assumption 1 1 . 1 1 . However, in the case when g = 0, the limit process 
obtained from the fluctuation central limit theorem will have variance that increases to 
oo as t — > oo. 

Finally, we show that in one particular case, the average of the limit process Z is driven 
by a Gaussian if-self-similar process R with H = > \. The driving process R is not 
fractional Brownian motion since it does not have stationary increments. This is directly 
related to the fact that the limit process Z satisfies Z(0) = and hence is not stationary. 
The process Z(T + •) is expected to become stationary as T — > oo. Similarly, the driving 
process R is expected to have stationary increments in the long run (i.e. R(T + •) — R(T) 
approaches a process with stationary increments, as T — > oo). We study the asymptotic 
behavior of the process Z(T + •) as T — > oo in Theorem 12.51 For simplicity we restrict 
here to the case b = q^_ 2 ) ■ It is shown that, as T — > oo, the process Z(T + ■) converges in 
distribution, in Cr[0, oo) to a stationary Ornstein-Uhlenbeck process driven by fractional 
Brownian motion with Hurst parameter H = > ^. 

The paper is organized as follows. We state all the results in Section Section [3] 
contains the proofs of Proposition 12.11 and Theorem 12.21 The proof of the central limit 
theorem will be provided in Section |U In Section [5] we represent the limit (centered) 
station average process Z as an integral with respect to a Gaussian process and give the 
proof of Theorem [23] on the moment stabilization property of the admission control policy 
g. Section 6 is devoted to the study of the asymptotic behavior, as T — > oo, of the process 
Z obtained from the central limit theorem, and the proof of Theorem 12.51 is given. 

The following notation will be used. We denote the set of non-negative integers by 
N and non-negative reals by R + . For a Polish space S, CgfO, oo) (resp. Dg[0, oo)) will 
denote the space of continuous (resp. RCLL) functions endowed with the local uniform 
(resp. Skorohod) topology. C will denote generic constants in (0, oo) whose value may 
change from one proof to next. 
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2 Model Formulation and Main Results 



We begin in this section with the evolution equations for the unsealed system. 
2.1 Unsealed System 

Let £o,ii i = I, ■ ■ ■ ,d, be independent Poisson random measures on [0, oo) x [0,oo) with 
common intensity rj = m x v, where m denotes the Lebesgue measure on [0, oo) and v is 
given in (jl.lip . Then, £o,i can be represented as 

oo 
3=1 

where < 5^1 < 5^2 < • • • are the jump times of independent unit rate Poisson processes 
for i = 1, . . . , d and Tjj are i.i.d. with distribution v. These Poisson random measures will 
be the building blocks for our counting processes iVj with desired intensities. 
Let / : R + x R — > R + be a function of the form 

f(t,y) = exp{-g(y-bt)}, (2.1) 

where g : R — > R is a function satisfying Assumption I l.li Let Xo = (-^0,1 > • • • ^o,d) T an d 
Yo = (^o,i; • • • > ^o,d) T be N rf and Rl-valued RCLL processes given through the following 
system of equations: 

00 

X 0>i (t) = N ,i(t) = & (i (So(t)) = HSi^Ao®}' 
Y ,i(t) = r A (t - 7d(s))6,i(d», dr) 

JB (t) 

= E ^'A(t-7o(^)), (2.2) 

j:S t:J <A (t) 

where 

A (t) = f f(s,Y (s))ds, Y (t) = JVy ,iW, 7o(t) = Aq 1 ^), (2-3) 

and 

B (t) = [0,Ao(t)] x [0,oo). 

Note that from Assumption 11.11 An is continuous and strictly increasing. Therefore, 70 is 
well defined and continuous as well. jo(Sij) is the jth activation time at the ith station, 
that is, the jth jump time of No t i(t). For jo(Sij) < t, t — jo(Sij) is the amount of time 
up to t since the jth session activation at the ith station and = 70 (S{ J) + Tij is the 
end time of the jth session at the ith station. Thus ry A (f — 70 (Sij)) is the work input 
by the jth activated source at the ith station, up to time t. 
From Assumption 1 1 . 1 1 it follows that 

f(t, y) = exp{-g(y - bt)} < exp{-g(-bt)}, V y > 0. (2.4) 
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In particular, / is a strictly positive function that is locally bounded, namely 

sup f(t,y) < oo, V T > 0. 
te[o,T],yew+ 

From this it follows that there is a unique solution to the system of equations (|2.2p - (j2.3p . 
Indeed, the solution can be constructed recursively between successive jump times of the 
underlying Poisson random measures {(,o,i, i = 1, • • • , d}. To see the basic construction, we 
consider for simplicity the case d = 1. The general multi-dimensional case can be treated 
similarly by arranging the jumps {Sij} in an increasing order. Simplifying notation, 
denote by < Si < S% < . . . the jump times of £ = £o,l and denote the corresponding 

session lengths by t\ , ti , 

For s G [0, Si), the solution to equations (12.2j) - (12.3p is given as 7o(s) = Aq 1 (s), where 
Aq(u) = J Q n f(v, 0)dv for all u > 0; and, writing iVo,i = N, Yb,i = Y, 

N(t) = Y(t) = for all t G [0,7o(Si)). 

Next, for s G [Si,S2), 7o(s) = Aj~ 1 (s), where 

Ai(«) = ( A ° (n) ' if UG [0 ' 7o(5l)) 
' \Si + / 7 "( 5l ) f(v,n A (t> - 7o(5i)))d«, if u > 70 (Si), 

and 

JV(t) = 1, y(f) = n A (t - 7 o(Si)) for all t G [ 7 o(Si), 7 o(S 2 )). 

One can now similarly write expressions for N(t) and Y(t) for i G [70 0%), 7o(Si+i)), for 
i > 1. 

This recursive construction shows in particular that 70 is a {J-~ u }-adapted process, 
where 

F u = a{ii{A) :AeB([0,u] x [0, 00)), i = 1, . . . , d}. 

Consequently, for any t > 0, Ao(i) = 7cT 1 (i) is a bounded {J^j-stopping time and there- 
fore, 

N 0>i (t) - A (i) = f 0> i([0,A (i)] x [0,oo)) 

is a {<5t} = {J 7 A (t)}-martingale, where £o,i = £o,« — is the compensated Poisson random 
measure associated with £o,i> i = 1, . . . , d. 

2.2 Scaled Workload and Main Results 

We now introduce the scaled system. For each fixed n G N, let £ n> \, . . . , £ n> d be independent 
Poisson random measures on [0, 00) x [0, 00) with common intensity measure 

r] n (ds,dr) = n a dsv n (dr), (2.5) 

where u n is introduced in (|1.12p . Define, for i = 1, . . . , d, 

n a n a 

Y n ,iit) = -^zt / r A(t- ln (s))^i(ds,dr), (2.6) 

n JB n {t) 



7 



where 



A«(t) = f f(s,Y n (s))ds, Y n (t) = ]-J2Y n>i (t), 7n (t) = A" 1 ^), (2.7) 
^ d i=l 

and 

B n (t) = [0,A n (i)] x [O,oo). 

Note that B n and 7n do not depend on i = 1, . . . , d. As for the unsealed system, we 
see that the solution (X n ,Y n ) T of the system (|2.6p exists and is unique on [0, oo) for 
each n, where X n = (X n>1 , . . . , X n4 ) T and Y n = (Y n>1 , . . . , Y n4 ) T ] and moreover, X n , 
Y n G D R d [0,oo). 

Consider the ODE 

ju(t) = af(t,U(t)), *>0, 



(2i 
C/(0) = 0, 



where 



/*CxD /"OO 

a = (/3 - 1)6* / r(6»r + l)^dr = (/3 - 1)6* / n 2 r(n9r + l)~^dr 
Jo Jo 



0(0-2) 

The following proposition will be proved in Section [3l 



Proposition 2.1 There is a unique continuous function U that solves (I2.8p . TTie solution 
satisfies 

sup 1 17 (t) - bt\ < oo. 

t>o 

In the case when b = a, we have U (t) = bt, for all t > 0. 

Remark 2.1 As an immediate consequence of the above proposition we have that f(t, U(t)) 
exp{—g(U(t) — bt)} is bounded above and bounded below away from 0, namely 

< M{f(t,U(t))} < sup{/(t, U(t))} < oo. 
*>° t>o 

Denote 

fy{t,y) = Q-f(t,y) =-exp{-g(y-bt)}g'(y-bt). 

Since g' is bounded from below and above, it follows from the above proposition that 
f y (t, U(t)) is also bounded below and bounded above away from 0, namely 

-oo < w£{f v (t, U(t))} < S up{f y (t, U(t))} < 0. 



Let X = (X 1 ,...,X d ) T = a- 1 (U,U,...,U) T , Y = (Y 1 ,...,Y d ) T = (U,U, . . . ,U) T . 
The following is the first main result of this work. 
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Theorem 2.2 As n — > oo, (X n , Y n ) T -> (X,Y) T in D R 2d[0, oo), in probability. Further- 
more, for any t > and any q £ [0, j3 — 2), 

sup n*|Y n (s) - f/(s)| 0, (2.9) 

0<s<t 

in probability, as n — > oo. 

Let V 6 Ck[0,oo) be given as the solution of 



V(t) = a f f y ( S ,U(s))V( S )ds-a6 2 -P f l^^ds. 



o jo 



(2.10) 



(2.11) 



From Remark 12. II the solution V of the above linear equation exists and is unique. 
Define 

z n ,i(t) = (r n 4t) - Yi(t) - ¥jj§£\ ,i = l,...,d. 

Our next result gives the limiting behavior of the processes Z n = (Z n> i, . . . , Z n 

,d) T - Note 

that Yi(t)+ V p l X is not the expectation of Y n ^ and hence, Z n> i in the above equation is not 
the conventional centered process of Y n j. However, from Proposition 12.11 and Theorem 
12.21 one can show Yi(t) = lim E(Y n j(i)). Also, as n increases to infinity, the term V S% 

n— >oo ' ' n 

tends to zero. Thus, the next result can be regarded as a central limit theorem for the 
scaled and (nearly) centered process Y n . 

Theorem 2.3 As n — > oo, Z n converges in distribution in D R d [0, oo) to Z = {Z\, . . . , Z^) T 
where Z satisfies 

Zi(t)= [ r A (t-*y(s))Ei(ds,dr) +a [ f y (s,U(s))Z(s)ds, i = 1, . . . ,d, (2.12) 



B(t) JO 

where 

rt 



A(t) = / f(s,U(s))ds, 7 (t) = A" 1 ^), (2.13) 



B(t) = [0,A(t)] x [0,oo), 
1 d 

z(t) = 1 j2 z i(t), 

and Ei,. . . ,£(2 are independent Gaussian random measures on [0,oo) x [0,oo) with com- 
mon control measure ds{(3 — l)6^~^r~^dr. 

Integrals with respect to Gaussian random measures characterized by a control mea- 
sure are defined, for example, in Chapter 3 of [TT] . 

Remark 2.2 When A = 1 (constant) or / = 1, note that A n (t) = t and 7 n (i) = t in 
(HTTP . B n {t) = [0,t] x [0,oo) and hence 

pt poo 

Y n,i(t) = — — r / / r A (i - (da, dr) 
n jo io 
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in (|2.6p . After the change of variables s — > s/n, r — > r/n, this can be written as 



1 



nt poo 

r A (nt — s)Cn t i (ds, dr ) , (2.14) 



where £ nj j is a Poisson random measure with intensity measure n a ~ - x ds {f3-l)9{9r+l)-Pdr. 
Written as (|2.14p . n a Y n ^ can be interpreted as the cumulative workload in the system 
scaled in time by n and where heavy-tailed workloads are associated with sources arriving 
at Poisson rate X n = n a . This is the view taken, for example, in |U|S]. It is well known 
that, after proper normalization and centering, the total workload converges to fractional 
Brownian motion in the so-called fast regime, that is, when 

n MT — > oo. 



This holds when a — (3 + 1 > 0, which is a part of our assumption (|1.15p . It is also known 
that the normalization of the right-hand side of (|2,14p (to the central limit theorem) is 

n 2 



(A^raM/ 3 " 1 )) 1 ^ n (a-/?+3)/2 
which coincides with that used in (|2.1ip . 

Remark 2.3 Let Z* = (Z\, ... , Z^) T be given as the solution of 



Z*(t) = R*(t) + a f f y (s,U(s))Z*(s)ds, (2.15) 
J o 

where 

R*(t) = / / (f(s, U(s))) 1/2 (r A (t - s))Ei(ds, dr), i = 1, . . . , d, 
Jo Jo 

and 

1 d 

z*(t) = -j2zm- 

One can check that R* = (R\, ... , R* d ) T and R = (R\, . . . , R d ) T have the same distribu- 
tion, where 



Ri{t) = r A (t - 7(a))E<(ds, dr), i = 1, . . . ,d. 

JB{t) 

Consequently, Z and Z* are equal in law and thus (|2.15p gives an alternative representation 
for the weak limit of Z n , as n — > oo. 

The following result shows the moment stabilization property of the admission control 
policy g. 
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Theorem 2.4 The following uniform moment bound holds: 

where /x := inf s >o{ — f y (s, U(s))} E (0,oo) andT(-) is the Gamma function. 

Remark 2.4 The case when there is no admission control corresponds to g = 0. Although 
the function g = does not satisfy Assumption ll.l| it can be shown along similar lines 
that in this case Theorem 12.21 holds with U(t) = at, and therefore, sup i>0 \ U(t) — bt\ will 
be finite if and only if b = a. Furthermore, Theorem 12.31 will hold as well (when b = a) 
but the moment stabilization property in Theorem 12.41 fails. 

Finally we consider the asymptotic behavior of Z(T + •) as T — > oo. Here we restrict 
ourselves to the case b = a. Then from Proposition 12. 1| (|2.12p and (j2. 13j) . the limit 
process in Theorem 12.31 can be written as 

ft poo rt 

Zi(t)= / rA(t-s)Ei(ds,dr)-K Z(s)ds, i = 1, . . . , d, (2.17) 



JO JO Jo 

where k = ag'(0) E (0,oo). 

Let Bh = {Bn{t), t > 0) be a standard fractional Brownian motion with Hurst 
parameter H = E (|, l) , namely, Bh is a mean zero Gaussian process with covariance 

E(B H (t)B H (s)) = i (t 2H + s 2H -\t- s\ 2H ) . (2.18) 
Let ^oo(O) be a normal random variable with mean zero and variance 

al— /9 />oo poo 

a 2 :=E(\Z oo (0)\ 2 )= d ^_ 2) J e- RV e~ KU \u - v\ 2 ~ p dudv < oo, (2.19) 

and let {Bh, Zoo(0)) be jointly Gaussian and the covariance function of Bh and ^oo(O) 
be 

al— (3 rt roc 

cov(B H (t),Z oo (0)) = ——— / e-™(u + v) 2 -Pdvdu, 
ad {(3 - 2) J J 



where a = J d((3-2)(3-i3)(4:-f3) • ^ e * ne fractional Ornstein-Uhlenbeck process given 

as the unique solution of 

Zoo(t) =Z oo (0)-k [ Zoo(s)ds + aB H (t). (2.20) 
Jo 

Theorem 2.5 Let b = a and let Z be as in Theorem \2.3l Then, as T — )• oo, Z(T + •) 

converges in distribution, in Cr[0, oo) to Z^ given by (|2.20p . Moreover, the process Z^ 
is stationary. 
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3 Law of Large Numbers 

In this section we will prove Proposition 12.11 and Theorem 12.21 



Proof of Proposition I2.lt Consider the ODE 

\u(t) = aexp{-g(u(t))} - b, t > 



(3-1) 

\u(0) = 0. V ; 

Clearly, a differentiable function u solves (|3.ip if and only if U(t) = u(t) + bt solves (|2.8p . 
From Assumption I l.H the function h(x) = aexp{— g(x)} — b, x E R, is locally Lipschitz. 
For each n G N, define h n (x) = h((x An) V (— n)), x£l, Since h n is a Lipschitz function, 
for any n E N, the ODE 

(u(t) = h n (u(t)),t>0, 

\u{0) =0 1 ' J 

has a unique solution u n . Let K be the unique solution of the equation 

aexp(-g(K))-b = 0, 

i.e. 5 (10 = log §. 

Then, for all n > |if |, if 6 > a, u n {t) < for all t and u n (t) decreases to K £ (— oo, 0); 
if 6 < a, u n (t) > for all t and u n (t) increases to K S (0, oo); and finally if b = a, 
u n (t) = for all t. Consequently, for any n > \K\, 

sup\u n (t)\ < \K\ (3.3) 

and u n solves (I3.ip . This proves the existence of solutions. 

Now consider uniqueness. Let u be another solution of (|3.ip . Let r = inf{t : > 
IX | + 1}. From unique solvability of (13.2p . for any n > |if| + 1, u(t) = u n (t) for all 
t G [0, r). From (j3.3|) we now see that r = oo. This proves unique solvability of ()3.ip and 
consequently that of (|2.8j) . Also, as noted above, 

sup |J7(i) - bi| = sup \u(t)\ < \K\ 
t>o t>o 

and U(t) — bt = u(t) = for all t, if b = a. The result follows. □ 

Next, we give the proof of Theorem 12.21 

Proof of Theorem 12. 2t Let £ n> i = £n,i — % be the compensated Poisson random measure 
associated with £ n> j, i = 1, . . . , d. Rewrite X n , Y n as 

X n ,m = \tnA B n(t)) + ^Vn(B n (t)) = — | -(^(t)) + A n (t), 
n « n a n a 

Yn,i{t) = — t I r A{t--y n (s))in,i(ds,dr) 

n JB n (t) 

f r A(t-jJs))n a ((3 - I)n9(n9r + l)~ p dsdr. (3.4) 

n ° JB n {t) 
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By the change of variables s = J f(u,Y n (u))du = A n (u), the second term on the right- 
hand side of (|3,4p equals 



t I'OO 



n 2 9{P-l) I I f(v,Y n (v)){r A {t - v)) {n9r + if 3 drdv. (3.5) 
Jo Jo 

Consider the inner integral in (|3.5p . For < v < t, by changing variables, we see that 
9(j3 - l)n 2 / (r A (t - a)) (n9r + l)~ p dr 



B{fi - l)n 2 
1 



o 

n(t— s) 



r(n6r + l)~ l3 dr + / (t - s)(nOr + l)'^dr 

Jt-S 

/•oo 

r(9r+l)~ l3 dr + n(t-s)j (9r + l)~^dr 

J n(t—s) 



1 



008 - 2) 0(/3-2)(n0(t-s) + iy»- 2 
Therefore, for each i = 1, . . . , d, 



a 1 



1 



(n9(t -s) + lf~ 2 



(3.6) 



Yn,i (^) 



1 



/r,a— 1 f 



+ / af(s,Y n (s))(l 



r A (t — 7n(s))| n ,i(rfs, dr) 
1 



(ra0(t - s) + l)/ 3 



-2 



and hence, 

X„,;(t) - Xi(t) 
Y nti {t)-Yi(t) 



rr 



—tn,i{Bn(t)) + / f(s,Y n (s))ds- / f(s,U(s))ds 



(3.7) 



(3.8) 



rnOt — 1 I 



r A (*- 7n(s))ln,i(^s,dr) 



+ / o[/(s,y n (s))-/(s,^(s))]d5-a 



o (n0(i - a) + If- 2 



ds. 



(3.9) 



We will first show (|2T9]h Since K n (t) > for all i > and /3 e (2,3), from J23|), we 
have that 



f{s,Y n (s)) 



-ds < 



1 



o (n9{t -s) + lf~ 2 ~ nP~ 2 9P- 2 J (t - sf~ 2 
Consequently, for every t > 0, 

f(s,Y n (s)) 



exp{-g(-bs)} < (sup < s < f expj-g(-feg)}) t 3 I 3 



n> 



'(3-/3) 



lim n q 



ds = 0, a.s., for all g £ [0,(3 — 2). 



n4oo'" y ( n 6)(t _ s ) + i)/3-2 
For n£ff, define the filtration {J~™} as 

7£ = <r{£;(A) : AGB([0,«] x [0,oc)),* = l,...,4. 



(3.10) 
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Then for each i = 1, . . . , d, £ n ,i([0, u] x [0, oo)) is an {J^j-martingale. As for the unsealed 
process in Section [2TTT 7„ is a continuous, strictly increasing {J^j-adapted process. Con- 
sequently, for every t > 0, A n (t) = 7^ x (i) is a {J^j-stopping time. Now consider the first 
term on the right-hand side of (|3.9p . that is, for i = 1, . . . , d, 



Observe that 

/(i) 



where, for i = 1, . . . , d, 



r A (t - Jn(s))£n,i(ds, dr) =: A n>i (t). 
A n , i (t)=ui 1 }(A n (t)), 



U { n}(u) = -—^ / r A (t - 7nO))+M<fe, dr). 

n i[0,u]x[0,oo) 



Note that U^ <(u) is a {J-"™}-martingale with predictable quadratic variation process 

2U.-D / {r A(t--f n (s)) + ) 2 i] n (ds,dr). 

™ 1 > V[o,«]x[o,oo) 

Using a change of variables, we have, for each t > and A n = (A nj i, . . . , A n ^) T , 
E(|A n (t)| 2 ) 



vi=l 



^ * V£l- / [°> A n(t)]><[0,°o) 



(r A (t — 7«(s))) 2 r/„(ds, c?r) 



< 



\i=l 

n9(/3 -l)d E / /■* Z" 00 ^ (s))(r A (f _ s))2( ^ r + ^ 
n V-/o jo 

n9 ^ ~J )d I I * exp{-g(-bs)}r 2 (n0r + l)-Pdrds 
n ° Jo Jo 



t poo 



+ n6(£_l)d i i exp |_ 5( _ 6s) |^_ s)2(n ^ r + 1) -/3 drds 



n JO Jt-s 



£ gg-ijd .u^ ^.{-rf-t.)}} |« /""-% r + 1) - Slfr(is 



JO 



0(/3 - l)d sup 0<s<t {exp{-g(-bs)}} 



rr 



(n(t- s)) 2 (8r + lY^dr ds. 

..(*—») 

(3.11) 
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For the first term on the right-hand side of (|3. 1 1 1) . note that 

ft rn(t—s) i ft rn(t—s) 



ft fn{i—s) i ft rn(t—s) 

- r 2 {9r + l)-^drds < — — / / r 2 ~^drds 

a Jo Jo 0Pn a J J 

n 3-(S-a f (j. _ S f-P ds 



n" Jo Jo V p n^ J J 

1 

F(3~-/3)' 

Since f3 — 2 < ^ + °~ 3 (or equivalently, a > j3 — 1), we obtain that, for all t > 0, 

1 2q ft fn(t—s) 



Iq rt pn(t—s) 

— / / r 2 (0r + l)-Pdrds^O, for all g G [0,/3-2). (3.12) 
n ° Jo Jo 



For the second term on the right-hand side of (|3.11j) . note that 



1 



t foo 

/ (n(t-s)) 2 (9r + l)-Pdrds 

Jn(t-s) 
t roo 



~ ek°l I, (n(t-s)) 2 r-Pdrds 

„3-P-a rt ^-P-a+A-P 

I (t-sf-Pds 



Thus we have 



n 2q ft f-oo 



n(t - s)) 2 {9r + l)~Pdrds -> 0, for all q £ [0, (3 — 2). (3.13) 

Jn(t~s) 



Combining (|3.1ip . (|3.12p and (|3.13|) we conclude that 

lim n 2q E (|A n (t)| 2 ) = 0, for all g € [0,0- 2). (3.14) 

n— >oo 

We argue next that n q A n = n q {A n ^\, . . . , A n ^) T converges to the zero process, in 
D R d[0,oo), in probability. In view of (|3.14p . it suffices to check that {n q A n } is tight. To 
prove tightness we will use a standard tightness criterion. Namely, we will show that for 
each fixed T > there exists Ct > such that for < h < 1 and h < t < T, 

n 4q E (\A n (t + h)- A n (t)| 2 |A n (t) - A n (t - h)\ 2 ) < C T h 2 . (3.15) 

The above inequality, together with the relative compactness of n q A n (t) for each t > 
(which follows from (|3.14|) ). yields tightness of {n q A n } (cf. Theorems 3.8.6 and 3.8.8 in 

&)■ 

Now fix T > 0. In order to show (I3.15D . it is sufficient to prove that, for any < h < 1 
and < t < T, 

n 4q E (\A n (t + h)- A n (t)\ 4 ) < C T h 2 . (3.16) 
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In the following, we use Ct > to denote a generic constant depending on T, 9 and f3 
whose value may vary from line to line. For r,s,h,t G R, denote 

<'*(r, s) = r A (t + h — 7n (s)) + - r A (t - 7 n(«))+- (3-17) 

Define, for i = 1, . . . , d, 

z4?(«) = ^zt / <'*(r, S )| nii (d a ,dr). (3.18) 

(2) 

Observe that ■ (u) is a { J^j-martingale with quadratic variation process 

n (a ' J[o,«]x[o,oo) v ' 
Since 7n (s) < i if and only if s < A n (t), we have 

A nyl (t + h)- A n>l (t) = U { *}(k n {t + h)). 

Recalling that A n (t+h) is a {^"j-stopping time, we have by the Burkholder-Davis-Gundy 
inequality that for some C > 0, 



E (\A n ,(t + h)- A n ,(t)\ 4 ) = E (|wg(A n (t + />))| 4 ) 




(3.19) 

Denote for r, s, h, t € R, 

#£'*(r, s)=rA(i + /i-s) + -rA(t - s)+. 
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By a change of variables and using (|A.ip from Lemma [A. II in Appendix we can bound the 
first term on the right-hand side of (|3,19p as 

2 I /' A i( t +' 1 ) r°° \ 
" E / / [$ h /{r,s)fn a e{l3- l)n{n9r + l)- p drds\ 



n 4(a-l) 

2 



(rt+h poo 
J J f(s,Y n (s))[^ t (r,s)]' i n a e( y f3-l)n(ner + l)-^drds 



n 



4(a-l) 

2 



o rt+h poo 

^ -^TTT / / e X p{- 5 (-6 S )}[^(r, S )] 4 n«0(^-l)n(n0r + l)-^rd S 



n 4 («-!) Jo Jo 



< C T rf- p ~ 6a h b - p . (3.20) 

For the second term on the right-hand side of f)3. 19|) , by a change of variables once more 
and using (1A.2P in Lemma lA.ll we have 



n4(a _ 1} -ixj j Q [$n\r, s)] 2 n a 9(f3 - l)n(n6r + l)^dr ds\ 



9 / / fA-n(t+h) fOO 

-E / 



2 



El/ / /(s,F n (s))[<' t (r,s)] 2 n a 0(/9-l)n(n0r + l)-Pdrds 



n4("-i) y \ ./„ Jo 

< 



2 S u Po < s <^{e-^(-^)} ^ _ + 

< Crn 6 - 2 / 3 - 2 ^ 2 ^). (3.21) 
Observing that 

P-2<- min{/3 + 3a - 5, 2/3 + 2a - 6} and min{6 - /3, 2(4 - /3)} > 2, 

and combining (|3.20p and (|3.2ip . we conclude that (|3.16p holds for every q £ [0, j3 — 2). 
This shows that for every such q, n q A n converges in probability to the zero process, in 
-D K d[0,oo) and thus for every t > 

sup n q \A n (s)\ —7-0, as n — > oo, in probability, for every q G [0, /3 — 2). (3.22) 

0<s<t 

Combining fl3~T0]) and (j3T22|) and recalling that \ Y? i= i{ Y n,i ~ Yi) = Y n -U, we have 
from that 

Y n (t) - U(t) = a f [f(s, Y n (s)) - f(s, U(s))] ds + S n (t), (3.23) 
Jo 



where 

f(s,Y n (s)) J -t m [* f(s,Y n (s)) 



o (n9(t - s) + I)/ 3 - 2 
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and n q S n converges in probability to the zero process, in Z) K d[0,oo). From Assumption 
we have that y \— > f(t,y) is a Lipschitz function on M + , uniformly in t on compact 
intervals, since 

sup \f v (t,y)\ = sup | - e W {-g(y - bt)}g'{y - bt)\ < Lexp{-g(-bt)} (3.24) 

J/GM+ 



for all t G [0,oo). Thus (I2.9P follows by an application of Gronwall's lemma to (I3.23p , 

Recall that for each i = l,...,d, £ n ^([0, it] x [0, 00)) is an {J^j-martingale and for 
every t > 0, A n (t) = 7 r ^ 1 (i) is a {J 7 " }-stopping time. Therefore, 



M£ht)=£ n , i ([Q,A n (t)] x [0,oo)) =£ n ,i(B n (t)) 



(3.25) 



is a {^™}-martingale, where Qf = F 1 ^ , t y By Doob's maximal inequality, for some C > 0, 



1 



— |M«(t)|>e) < 



o<s<t n L 



CE|Mff(t)| 2 

n 2a e 2 

CE(rj n ([0,A n (t)] x [0,oo))) 



< 



CE (ftf(s,Y n ( S ))ds 
n a e 2 

c fo ex P{-9(-bs)}ds 



n a e 2 



Combining (|3.26[) and (|3.25p we have that, as n — >• 00, 



P sup 



o<s<t n 



in,i(B n {t)) 



> e 



0. 



(3.26) 



(3.27) 



Thus the first term on the right-hand side of (|3.8p converges to the zero process, uniformly 
on compacts, in probability, as n — > 00. Finally, combining (|3.24p . (|3,27p . (|3.10p and (|3.22l) 
(with q = 0), we deduce applying Gronwall's lemma to (|3.8p and (|3.9p that (X n , Y n ) T — > 
(X, Y) T in D R 2d[0, 00), in probability, as n — > 00. □ 



4 Central Limit Theorem 

In this section we prove Theorem 12. 3[ From (|2.10p and (|3.9|) . we can write (|2.11|) as 



Z n ,i{t) 



B n (t) 



(r A (t - 7„(s)))S n)i (ds, dr) 



p-2 



ds 



+a / n 3 *^ f(s,Yn(s)) - f(s,U(s)) - f y (s,U(^ {s) 

a +.fl-3 /* af(s,Y n (s))ds c.-0+i f* a0 2 -Pf(a,U(8)) ds 
~ U 2 J (n9{t-s) + iy- 2+n 2 J {t-sf~ 2 ' 



(4.1) 
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where, with a n = n^ a ^ +1 )/ 2 7 



y+P-3 



^n,i(A) — n a—l &m("^-) ~~ °ri £rM(^) 

is a random signed measure on [0, oo) x [0, oo), i = 1, . . . , d. Note that 

Var(E n) i(A)) = vP^m x u n {A), i = l,...,d, for A e B(R+) with m x f n (A) < oo. 
Note also that 

f{8, %{S)) - f( S , U( S )) - f y {8, U(S))^ 

= (Y n (s)-U(s)) f [f y (s,U(s) + x(Y n (s)-U(s)))-f y (s,U(s))}dx 
Jo 

+ (Y n (s) - U(s) - ^g) /„(,, U(s)). (4.2) 
Thus the middle term on the right-hand side of (|4.1|) equals 

a n ^{Y n (s)-U{s)) [f y (s,U(s) + x(Y n (s) -U(s))) - f y {s,U{s))]dxds 



JO 



+ a f y (s,U(s))Z n (s)ds, 
Jo 

where, recall, Z n (s) = \ Ya=i z n,i{s). Let 



Rn,i( t )= (r A (t — 7 n (s)))£ n) j((fe, dr), i = 1, . . . ,d, (4.3) 

JB n (t) 



ft fi 
C n (t) = a n B± ^(Y n {s)-U(s)) / [f y (s,U(s) + x(Y n (s) - U(s))) - f y (s,U(s))]dx ds, 
Jo Jo 

(4.4) 

Letting 7^ n (t) = (fc n ,i(t), . . . , fc n ,d(t)) T , where 7^(4) = R n ,i{t) + C n (t) + -D n (t), we can 
rewrite equation (|4.ip as 

Zn,i(*) = 72n,i(t) + a [ f y (s, U(s))Z n (s)ds, i = l,...,d. (4.6) 
J 



Proof of Theorem 12.3b Define ^ : D K d[0, oo) — >■ D R d[0, oo) by 

[V>(x)]»(i) = a?i(*) + a / fy{s,U(s))i()(x)(s)ds, i = 1, . . . , d, x € D U d[0, oo) 
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where ip(x) = 3 ^f=i[V'( a; )]i- Then, V is a continuous mapping from D R d[0,oo) to 
D R< {[0, 00). Also, from (|4.6j) we see that Z n = {Z n ,\\ ■ ■ ■ 1 Zn,d) T = ipffin)- 

Combining Lemmas 14 . 1 114 . 2 1 and 14 . 3 I below . we see that TZ n converges to R = (R\, . . . , Rd) T , 
in distribution, in D R d[0, 00), where 



Ri(t)= (r A (t — j(s)))Tii(ds, dr), i = 1, . . . ,d, (4.7) 

JB(t) 

and Sj, i = 1, . . . , d, is as in Theorem 12.31 The result now follows from continuous map- 
ping theorem. □ 

The next three lemmas were used in the proof of Theorem 12.31 above. 

Lemma 4.1 Let C n be as given in (|4,4p . As n — > 00, sup 0<s<i |C n (s)| — > 0, in probability, 
for every t > 0. 

Proof: From Assumption 11.11 we have 

\fyy(t,y)\ =| exp{-g(y - bt)}[g'(y - bt)} 2 - exp{-g{y - bt)}g"{y - bt)\ 
<{L 2 + L)e- g{ y-^ < (L 2 + L)e- 9{ - bt) =: c(t) 

for all y £ [0, 00). Consequently y >— > f y (t,y) is a Lipschitz function on R + , uniformly in 
t in compact intervals. Therefore 

\C n (t)\ <a J c(s) (n^^\Y n {s) - U{s)\) 2 ds. (4.8) 

The result now follows by noting that Q+ f~ 3 < (3 - 2 (see (fl~T5]) N | and using $ZJft. □ 



Lemma 4.2 Let D n be as given in (|4.5p . As n — > 00, supo< s <t |.D n (s)| — > 0, in probabil- 
ity, for every t > 0. 

Proof: Note that 

DJt) _ ^(j;*^^-^ ^*^ * 

+ n 2 ^ /" , ° / "' Cf(S) ' - n"- 2 /' , y { '' %{s) l .,1- I . 







(n0(t - s) + If- 2 J (nd(t - s) + If- 2 

(4.9) 
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For the first term, note that 

< ?a&-l*f{fiM*)) 



< 



< 



(t - sf- 2 
a6 2 ~Pf( S ,U(s)) 

a6 2 -? 



ds - n 13 ' 2 
( 1 



af(s,U(s)) 
o (n9(t - s) + lf~ 2 

1 \ 



ds 
ds 



sup {exp{- 5 (-6 S )}} I t 3 ^ - (t + — a f-P + ( ^ 



3 — P o - < 



\3-/3 



1 \ 



3-/3^ 



3-/3 <^ {eXP{ - 5( - 5s)} H^ 



3-/3 



Also, from the Lipschitz property of / (see (|3.24p ) we have 
< l ^— -^-ds 



n 



o {n9{t-s) + iy 



< 



P- 2 , . 

o (n9(t - s) + If- 2 

t a6 2 -PLe-^- bs ^\Y n {s) - U(s)\ 

o (t-s + ±y- 2 

a0 2-f3 Lt 3-f3 

sup {exp{-g(-bs)}} sup {\Y n (s) - U(s)\}. 

-J — P 0<s<t 0<s<t 



-2 



ds 



Combining (|4.9[) - f)4. lip , we have 



" Q+S-5 a— 0+1 

sup \D n (s)\<C + n^^\U n {s)-U{s) 

0<s<t L 



(4.10) 



(4.11) 



(4.12) 



From (fTT5D we see that SE±|z§ < o and 2=|±i < /3 - 2. The result follows using □ 



Lemma 4.3 Let R„ = (R n ,i, • • • , R n ,d) T and R = (^i> • • • > -^d) T ^ e gu>en fry (|4.3|) and 

(|4.7p respectively. As n —> oo, R n converges to R in distribution, in D K d[0, oo). 

Proof: Let R n = (i? n ,i, • • • , Rn,d) T , where 



i2n,i(t) = / [r A (t - 7(s))]E n , i (ds J dr), t > 0. 
Since A n (t) is an {J^j-stopping time for each t > 0, l[o,A„(t)](s)[r A (i — 7n(s))] is J 7 " 
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predictable. Thus applying the isometry property of the stochastic integral, we obtain 

2 



E (R^i® - R n .i(t) 

= E ( / (l [0iAn(t)] (s)[r A (t- 7n (s))] - l [0 ,A W] (s)[r A (t — 7 (s))]) E n ,<((fo,dr) 

JO JO 

(/»oo /»oo \ 
^ ^ l [0 ,A nW ](s)[rA(t- 7n ( S ))-rA(t- 7 ( S ))]V- 1 i/ ri (dr)d S J 

/>oo 

+2E(|A n (t)-A(t)|) / (rAi)V~V(dr) 



= 2E( / jf i [0Ai(t)] ( s )[ rA (t_ 7n ( s ))_ rA (t- 7 ( s ))] 2 ^-^dr^ 

/■oo 

+2E(|A„(i) - A(t)|) / (r At) 2 6((3 -l) n P(n9r + l)-Pdr. (4.13) 
J o 

By the dominated convergence theorem, we have by using the fact that f3 £ (2,3), 

/■00 /'OO 

lim / (r f\t) 2 9{fi — l)nP (n8r + dr = \ (r A t) 2 9 l ^^(j3 — l)r _/3 <ir < oo. (4.14) 
n ^°° Jo 7o 

Note that 

< min{A(t), A n (t)} < max{A(i), A n (t)} < / exp{-g(-bs)}ds. 

Jo 

Consequently 

|A n (t) - A(t)| < 2 / exp{-g(-bs)}ds. 
Jo 

Also, from (|3.24|) we have 

\A n (t) - A(i)| < L sup |tf n (s)-i7(«)| / exp{-g(-bs)}ds. 

0<s<t Jo 

Thus (|2.9h and the dominated convergence theorem yield 

lim E(|A n (t) - A(t)|) = 0. (4.15) 

n— too 

Combining (|4.14|) and (|4.15|) we have that the second term on the right-hand side of (|4.13|) 
converges to as n — > oo. 

Now we consider the first term. From the definitions of 7n and 7 we see, for any s > 0, 

7n(s) _ ff(s) 

f(z,Y n (z)dz= / f(z,U(z))dz. 
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Consequently, 

n( s ) 



7n(s) 



f(z,U(z))dz- / f(z,U(z))dz 



7n(s) 



7n0) 



f(z,Y n (z))dz- I f(z,U(z))dz. 
Jo 

(4.16) 

Since f(z,U(z)) is bounded below away from (see Remark 12. ip . there exists a c > 
such that 



7 (s) 



f(z,U(z))dz 



7n(s) 



f(z,U(z))dz 



>c|7(a)-7n(*)|. (4-17) 



On the other hand, from (|3.24p . we obtain that for any s < A n (t) (equivalently, 7 n (s) < £) 



7n(s) 



/(z,y n (z))dz 



7n(s) 



f(z,U(z))dz 

7n0) 



< L sup |y n (2;) — ^7(^)1 / exp{— g(— bu)}du 

0<z<7„(s) 

< L sup \Y n (s) — U(s)\ I ex.p{—g(—bu)}du. 

0<s<t 

Combining (|4.16p - (|4.18p we have that 

1[0,A„(*)]( s )It«( s ) - l(s)\ < — sup \Y n (s)-U(s)\ I exp{-g(-bu)}du. 

C 0<s<t 



Using (|2,9p we now obtain 



l[O,A n (t)](s)|7n(s)-700! -> 



(4.18) 



(4.19) 



(4.20) 



in probability, as n — > oo. An application of the dominated convergence theorem now 
shows that 



lim E 

n— >oo 



' l[o,A„W]( s )[ r A (* " 7nW) - r A (t - 7(«))] 7 

V 



0(/3 - l)r/ 
nOr + l)/ 3 



<ir ds ) = 



Thus the first term on the right-hand side of (|4.13p converges to as well. 

Combining the above observations we have that for each i = 1, . . . , d and t > 



(4.21) 



lim E [Rn,i(t) — R n .i(t) ) =0. 



(4.22) 



Note that 



Rn.i(t) 



n 2 1 l B( t)(s,r)(r A(t->y(s)))£ n ,i(ds,dr). 



'[0,oo)x[0,oo) 

For each fixed i = 1, . . . , d, let us show the weak convergence of the finite-dimensional 
distribution of R n ^. For any < t\ < ■ ■ ■ < tk < oo, denote f n (s, r) = (/™(s, r), . . . , f£(s, r)) T 



where / n (s,r) = n 2 l B / t n(s, r)(r A (tj — j(s)), j = 1, . . . ,k. Then 



[0,oo) x [0,oo) 



f^(s,r)£ nti (ds,dr). 
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lim / 



One can show by the change of variables that 

fJ(s,r)fl l (s,r)7 ln (d S ,dr) 

[0,oo) x [0,oo) 

lim f f^(s,r)fl l (s,r)n a (/3 - I)n9(n9r + l)"? dr ds 

n ^°°J [0,00) x [0,00) 
rtjhti poo 

lim / / f(s,U(s))[r A (t,- - s)][r A (tj - s)]n /3 (/3 - lWnflr + l) - ^drds 
n ^°°Jo Jo 

ftjAti poo 

f(s, U(s))[r A (tj - a)] [r A (i, - a)] (/3 - l)0 1_ ^r - ^ 'dr ds 
= E(Ri(tj)Ri(ti)). (4.23) 

P-a-l g— q— 1 

Since |/ n | < n 2 £ fc an d lim n 2 = 0, we deduce that lf\fn\ >£ \ = for large enough 

n— too 11 1 ; 

n, and hence, for each e > and j, if n is large enough, 

1 {|/™|> £ }l/j 1 (' s ' r )| 2r ?™( d ' s ' dr ) = °' 

[0,oo) x [0,oo) 

From Theorem 6.1 in [6] it now follows that (Rn^ti), . . . , R n ,i{tk)) T =^ {Ri{ti), • • • , Ri{tk)) T 
as n increases to infinity, for each i = 1, . . . , d. Since R n has independent components, we 
have that the finite-dimensional distributions of R n converge to those of R. Using (j4.22j) , 
we then obtain that the finite-dimensional distributions of R n converge to those of R. 

It thus suffices to show that {R n } is tight in D R d[0, 00), for which, it suffices to prove 
the following estimate: for each fixed T > there exists a constant Ct > such that for 
0</i<landO<t<T 

E (|R n (t + h)- R„(t)| 4 ) < C T h 2 , (4.24) 
Recall the definition ofu£} in (pTTHjl . Then 

RnAt + h )~ Rn,i(t) = n2± ^Un}( A n(t + K)) = {A n>i {t + h) - A n ^{t)) . 

From (|515|> . (|3T20l) and (|3T2ll . we now have 

E(|R n (t + /i)-R n (t)| 4 ) = n ii ^ 1 E(\A n (t + h)- A n (t)| 4 ) 

< c - n 2a+2 l 8-6 ^5-/3-3^6-/3 + n 6-2/3-2a ^2(4-/3)^ 

= C (rr^ +1 ^ 6 ^ + Z^ 4 "^) < C/* 2 , (4.25) 

where the last inequality follows from a > (3 — 1 and 2 < (3 < 3. This proves the desired 
tightness and the result follows. □ 
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5 The Moment Stabilization Property 

In this section, we will prove Theorem 12.41 Let Z be as in Theorem 12.31 and let R = 
(R\, . . . , Rd) T be the Gaussian process introduced in (|4.7p ■ Then Z = ^ Yli=i %i satisfies 



Z(t) = R(t) + a f y (s,U(s))Z(s)ds 



(5.1) 



where R = \ Y^i=i^i- Note that R is a zero mean Gaussian process. We begin by 
computing the covariance functions of R4, i = 1, . . . , d, and R. 

Lemma 5.1 The covariance functions of the Gaussian processes Ri, i = 1, . . . , d, and R 
are given respectively by 



cov{R l {s),R l {t)) = E(Ri(s)Ri(t)) 



and 



cov(R(s),R(t)) = K(R(s)R( y t)) 



,1-/3 



s ft i-uAv 

Jo Jo 



exp{—g(U {z)—bz)}{u\/v—z) 1 13 dz du dv, 

(5.2) 



s ft ruAv 



JO Jo 



for any s, t > 0. 



exp{— g(U(z)— bz)}(u\/v— z) 1 dzdudv, 

(5.3) 



Proof: Without loss of generality, we assume that < s < t < 00. For each i = 1, . . . , d, 
by the change of variables u = Jq exp{— g(U(x) — bx)}dx, we have 

9^- 1 E(R i (s)R i (t)) 

[r A (s - j(u))][r A (t - j(u))]((3 - \)r~^drdu 



B(s) 

S /"OO 



J 



exp{-g(U(z) - bz)}[r A (s - z)][r A (t - z)]{(3 - \)r~^drdz 



Jo 



+ 



Js-z 
s roo 



exp{-g(U(z) - bz)}r 2 (/3 - l)r~ p drdz 

S ft — Z 

exp{— g(U(z) — bz)}(s — z)r(/3 — l)r~^dr dz 
exp{-g(U(z) - bz)}(s - z)(t - z)(/3 - l)r~ p drdz 

Jt-z 

exp{-g(U(z) - bz)} 

^(s ~ z?-P + JL( a - z)[(s - zf-P - (t - zf-?\ 



dz. 



(5.4) 



A simple calculation shows that 
3^(* - z)*-P + ^ (s - z) \{s - z?-P -(t- zf-P] 



s rt 



(uVv-z) 1 p dudv. (5.5) 
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Combining (|5.4p and (|5.5p and by changing the order of integration, we obtain 
fl^E^s)^)) = / exp{-g(U{z) -bz)} (uV v - z^dudv dz 

JO J z Jz 

/•s rt rul\v 

= exp{-g(U(z) -bz)}(uV v - z) x ~^dzdudv. 

Jo Jo Jo 

This proves (|5.2p . Equation (|5.3p is now immediate on noting that J?j, . . . , i?^ are i.i.d. □ 

In the next lemma, we give a bound on the second moments of the increment of the 
Gaussian processes R4, i = 1, . . . , d, and R. 

Lemma 5.2 For any s,t > 0, the following bounds hold: 
and 

»(i*w-A(.)i')< j(< ,_^ 1 l^ (4 _ ffl (t-.r g . (5.7) 

where K\ := sup s>0 {exp{— g(U(s) — bs)}}. 

Consequently, Gaussian processes R\, . . . , i?<i, -R frove versions that are Holder contin- 
uous of any order p e (0, (4 — (3)/2), on [0,T], for all T > 0. 

Proof: Fix < s < t < 00. From Lemma 15.14 for each i = 1, . . . , d, 
E(\Ri{t) - Ri(s)\ 2 ) = cov(Ri(t), Ri(t)) - 2cov(Ri(s), Ri(t)) + cov(Ri(s), Ri{s)) 

= e 1 -? / / / exp{-g(U(z) -bz)}{uyv-zf- p dzdudv 

J s Js Jo 

/t pt fuAv 
/ / (uVv-z^-^dzdudv 
Js Jo 



Ki9 



1-/3 r-t r t 



0-2 



du dv 



Kie i-0 rt rt 



< — — / I \a — c\" ^dudv 



[u-v 



P-2 

2K^6 l 



(/3-2)(3-/3)(4-/3) 



2 ~ p dudv 



(t-s) 4 -?. (5. 



This completes the proof of (|5,6p . Inequality in (|5.7p is now immediate. The second 
statement in the lemma now follows from Kolmogorov's continuity criterion. □ 



Proof of Theorem 12.41 relies on an explicit representation for the solution of equation 
(|5.ip . For that we begin with an indefinite integral of a deterministic function with respect 
to the Gaussian process R. 
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Denote by £ the linear span of indicator functions of the form l( s ,t] '■ ^+ ~~ * ^> 
< s < t < oo. Consider the inner product on £ given by 

(l {0 , s] ,l m )n n = cov(R(s),R(t)). (5.9) 

We denote by 7-Lr the Hilbert space obtained as the closure of £ with respect to this inner 
product. Define R : £ -> L 2 (0, F, P) as 

H(l (0)t ])) = R(t), < t < oo, 

where the definition is extended to all of £ by linearity. Clearly E(|.R(</>)| 2 ) = (cp, faiUn f° r 
all (j) € £ . We can now extend the definition of R to all of %r by isometry. Occasionally, 
we will use the notation 



[■OO 

R{fa) = / <j>(t)dR(t), 4> 6 ftfi. 

JO 



Define 

p(u> „) = J ^ eM-9(U(z) - bz)}(u V v - zY-Pdz, if « / t,; (g 1Q) 

1 0, if u = u, 

for all u, w > 0. From Lemma |5.1[ we see that 



cov(R(s),R(t)) = [ [ p{u,v)dudv (5.11) 

J ■/ 

and for any cj), cf> £ £ 

(Mu s = / (f>(u)4>(v)p{u,v)dudv. (5.12) 
Jo Jo 

It can be shown that 7~Lr contains all measurable functions <j) on M + satisfying 

\4>(u)\\<j)(v)\p(u, v)dudv < oo. (5.13) 



One can also show that equality fj5. 12f) holds for 4>, <j> that satisfy (I5.13p . 

This type of isometry is considered in |9] (see Chapter 5) and [1] with respect to 
fractional Brownian motion and general Gaussian processes respectively. 

Remark 5.1 If <p : [0, oo) — > R is continuous, then, for any t > 0, the function fa defined 
by fa(') = l[o,t](0^(") satisfies (I5.13p . Consequently, fa is in 1-Lr and we write formally 

R{fa){t) := R(fa) = f cj)(s)dR(s). (5.14) 
J o 

Remark 5.2 If </>(•) is Holder continuous of order pi > 1 — on [0, t], for every i > 
0, as a result of Young integration theory |12j . the pathwise Riemann-Stieltjes integral 
Jq 4>(s)dR(s) exists, since R is Holder continuous of any order p £ (0, ^5^). Zahle [13] 
showed (see Proposition 4.4.1 therein) that R((f>)(t) is Holder continuous of the same order 
as R on [0,T], for every T > 0. The indefinite integral R(cf))(-) on the right-hand side of 
(|5.14p coincides with the pathwise Riemann-Stieltjes integral. 
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We now proceed to the proof of Theorem 12.41 



Proof of Theorem 12.4b Define 

4>(t) :=expj-a^ f y (z,U(z))dz^,4>(t) = explaJ f y (z,U(z))dz 

Then the derivatives of <j> and <j) are 

<f/{t) = -af y (t,U(t)Mt), 4>'(t) = af y (t,U(t))4>(t). 

Remark 12,11 implies that (/)' and <j)' are bounded on any compact interval, and hence, (f> 
and 4> are locally Lipschitz continuous. From Remark 1 5. 2 \ the indefinite integral 

H(0)(i) = jf cj)(s)dR(s) = J expj-a^ f y (z,U(z))dz\ dR(s) 

is well defined as a Riemann-Stieltjes integral, and for every T > 0, R(<p)(t) is Holder 
continuous on [0, T] of any order p G (0, (4 — /3)/2). 
It follows from Theorems 3.1 and 4.4.2 in [131 that 



4>(t)R{<p){t) = f 4>(s)(f>(s)dR(s)+a I f y (s,U(s))4>(s)R(<f>Xs)ds 
Jo Jo 

= R(t) + a[ f y (s,U(s))fa)R((f>)(s)ds, (5.15) 
Jo 

which implies that 4>(t)R((J))(t) solves (|5.1|) . Thus, the solution Z to (|5.1|) can be written 
explicitly as 

Z(t) = <j){t)R{(j)){t) =exp|a J f y (z,U(z))dz\ expj-a J f y (z,U(z))dz\ dR(s). 

(5.16) 

By the isometry of the mapping R we have, on letting (j>t(u) = 0(u)l[ot](w), 

f'OO f'OO 

E(|Z(t)| 2 ) = \4>(t)\ 2 / / <t> t {u)<f>t(v)p(u,v)dudv 
Jo Jo 

= \4>(t)\ 2 / / <f)(u)<f)(v)p(u,v)dudv 
Jo Jo 

= J J exp|a J / B (*,ff(z))(fe}exp{a J f y (z,U(z))dz^j p(u,v)dudv. 

(5.17) 



2S 



Recall the definition of \x from the statement of Theorem 12.41 From Remark 12. 1\ fj, £ 
(0,oo). Then by a calculation similar to (|5.8p . we have 



t r t 



Jo 



e^EdZit)] 2 ) < e -a^t-u) e -a^(t-v) p (u, V )dudv 



t rt ruAv 

/ / e-^-^e-^-^iuV v - z) l -Pdzdudv 
Jo Jo 



< un 2 . f f e-^-^e-^-^iu-vf-^dudv 
~ d(/3 - 2) J J v V 1 



t r t 



< 2 / / e - a K*-«)( n _ u )2-/3 dndt; 



d(/3 - 2) 7 7, 
2 

d(/5-2)(3 - /3) 7 



e-°"f*- , ')(t-«) 3 -^du 



S d(/3-2)(3%) (a< .)^ r(4 - ffl ' (5 ' 18) 



for all f > 0. □ 



6 Fractional Ornstein-Uhlenbeck Process 

We now proceed to the proof of Theorem 12.51 in this section. Throughout this section we 
take b = a. From Proposition 12.11 it follows that 

U(t) = bt, f(t, U(t)) = 1, f y (t, U(t)) = -g'(0), for all t > 0. (6.1) 

For notational simplicity, we will only present the proof for the case 8 = 1. 
In this special case, the SDE f)5. 1 j) can be written as 

Z(t) = R(t) -at g'{0)Z(s)ds = R(t) - « / Z(s)ds, (6.2) 
Jo Jo 

where k = afj, = ag'(0) > 0, and 



1 P* poo 

R V = lH / / rA(t-«)S 
d i=l Jo Jo 



(ds, dr), 



for any t > 0. 

From Lemma |5. 11 we have that the covariance of -R is given by 



1 



cov{R(s),R(t)) =E(R(s)R(t)) = -111 (uV v - z^'^dz dudv 



s ft ruAv 



JO JO 
1 



d(J3-2)J J 



s rt 



(\u- v\ 2 ~P - (uVv) 2 -P)dudv, 

(6.3) 
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and from Lemma 15. 2} we recall that the sample paths of the process R are Holder contin- 
uous on [0, T] of order p, for any p £ (0, ^^). 

Recalling the definition of the indefinite integrals with respect to the Gaussian process 
R, the solution of the SDE in (16. 2\\ can be explicitly written as 



r-t 

Z(t) = e~ Kt I e KS dR{s). (6.4) 
J o 

We now consider the asymptotic behavior of the process Z(t), as t — > oo. For T, t > 0, 
let fl T (i) = R(T + t) - R(T). From ([62]), we can write 

Z(T + t) = Z(T) + R T (t)-K I Z(T + s)ds. (6.5) 

JO 

Recall the parameters <Tq and a introduced before Theorem 12.51 (also recall that 6 = 1). 
Lemma 6.1 We have 

lim E(|Z(T)| 2 )=a 2 . (6.6) 

T->oo 

Proof: From (I6.3p . (|6.4p and the isometry property of of the mapping R, we have 

p -2kT I 

d{(3-2) 

From (f5TT8D . 



E(|Z(T)| 2 ) = - T7 J—^ e - 2KT I I e KV e KU \u - v^ 13 - (uVv) 2 '^ du dv. (6.7) 

11 Jo Jo 



for all T > 0. Also, by a change of variables, we see that 

h x (T) = f f e- KV e- KU \u-v\ 2 -Pdudv. (6.9) 

d (P - 2) 7 Jo 

Since the integrand is non-negative, the function h\(T) is increasing in T. Together with 
the boundedness of h\, the limit lim h\(T) exists and is finite. Furthermore, 

T^oo 
-1 /'OO /*oo 

r^Pl-S^y. /„ «-«->-•!*-''*.* = <* («0) 
Next, for the second term on the right-hand side of (|6.7|) . note that 



1 



T ,-T 



h 2 (T) := d{p _ 2 f ~ 2KT I I e KV e KU (uVvf-Pdudv 
2 



d(P-2) 
2 



o Jo 

c- 2hT I ' f V e KV e KU v 2 -Pdudv 
o Jo 

■T 



dK((3-2) J 



30 



An application of l'Hopital's rule shows that 

lim h 2 (T) = — — — - lim l - = 0. (6.12) 

Using (16. 1QH and (I6.12p . we now have that 

i poo roo 

lim E(\Z(T)\ 2 ) = — -/ / e- KV e- KU \u-v\ 2 -Pdudv = al (6.13) 

T^oo d (p - 2) J J 

The result follows. □ 

In the next lemma, we compute the limits of second moments and the covariances of 
Rt(-), as T — > oo. 



Lemma 6.2 For any t > 0, 



.Urn E(\R T (t)Y) - d(tj _ 2)(3 _ m _ 0) = oV, (6.14) 



and, for any t > s > 0, 



2 

lim E(R T (t)R T (s)) = t^l 3 + s 4 "^ - (t - s) 4 "^ . (6.15) 
Proof: For any t, T > 0, note that 

i? T (i) = R(T + t) - R(t) = I + l {T>T+t] (s)dR(s). 

Jo 

By the isometry property of the mapping R and a change of variables, we have 

lim E(\R T (t)\ 2 ) 
T-»oo 

I rT+t rT+t 

= lim—- -/ / I \u- v\ 2 ~P - (u Vd) M dudv 

T^oo d (J3 - 2) J T J T V 1 v 1 J 

T+t 

lim , .," / / ((« -u) 2 "' 3 -v 2 ^] (iucfo 

r->ood(/3-2)(3-/3) 7 T V 7 \ yj j 

1 (v 3 ~P - (3 - /3)(y + T) 2 ~Pt) dv 



t^oo d ((3 — 2) y T y T 

2 



lim 

T^ood(/3-2)(3-/3) 7 



2 



d(/3-2)(3-/3) y 



d(/3-2)(3-/3)(4-/3) 
This proves (|6.14|) . 



v 3 13 dv 
— - a 2 t A -P. 
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For any t > s > 0, by the isometry property of the mapping R and a change of 
variables, we have 

lim E{R T (t)R T (s)) = lim E((R(T + t) - R(T))(R(T + s) - R(T)) 

T->oo T->oo 



1 



lim , , 

T^oo d(fi -2) J T J T 

-y r T+s r T+s 

lim 



u — V 



mVu 



,2-/3 



du dv 



T->oo d(f3 — 2) J T J T 



u — v 



+ lim 



T->oo d((3-2) J T+S J T 



T+t r T+s 



2-/3 - (iiV«P) ciucfo 
{{v - u) 2 - p - v 2 -^ dudv 



lim E(\R T {s)\' 

T->oo 



+ lim 



1 



r^ood(/3-2)(3-/3) 



/ T+ * ((v - Tf-P — (v — T — sf~P - (3 - /3)w 2 - /3 s) cfo 



+ 



d(/3-2)(3-/?)(4-/3) d(/5-2)(3-/3) 



1 



d(/3-2)(3-/3)(4-/3) 

2 



(t-s) 



4-/3 



~2 



s 4-4» +t 4-/9_ (t _ a) 4-0 



where the fourth equality follows from (|6.14|) and a change of variables. This proves (|6,15[) 
and the result follows. □ 



The next lemma gives the limit of the covariance function of Z(T) and Rxit), as 
T ->• oo. 

Lemma 6.3 For any t > 0, 

lim aw(Z(r), £t(*)) = lim E(Z(T)i2 r (t)) = / / e"^(n + t,) 2 "^ d«. 

T^oo T^oo d{p-2) J J 

(6.16) 



T+t 



l(T,T+t](s)dfl(«) 



Proof: For any i, T > 0, note that 

fl T (t) = fl(T + i) - i?(T) 



Then, it follows from the isometry property of the mapping R and (|6.4p that 
lim cov(Z(T),R T (t)) = lim E(Z(T)fl T (i)) 

T— >oo T->oo 



rid(F2) y / T 

lim 



((«-«) 



2-/3 _ „2-j8 



u p ) dudv 



T^oo d{p-2) 



JT 



'(«-t;) 2 "^d«dt> 



lim , , 

T->oo d(/3 — 2) J J T 



e KV u z - p dudv. 



(6.17) 
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For the second term on the right-hand side of (|6.17p , we have 

— kT pT rT-\-t j. —kT rT 

lim 6 / / e KV u 2 ^dudv < lim 6 - / e KV T 2 ^dv 

d{fi-2) J J T t^oo d(p-2) J 

+(l - e~ KT ) 
= lim —4- — = 0. (6.18) 

For the first term on the right-hand side of (|6.17|) . using a change of variables, we obtain 

-kT rT r T+t 

ttod(p=2)J .//' 

-kT rT 

4 im 77r ovi m / e ^ T + f ~ v ? ~ ( T " v ^ dv 
T^ood{f3 - 2)(3 - P) Jo 



lim — - / / e KV (u-v?- p dudv 



lim "777^ w ttt / e^tf* + l-) 3 ^ - U 3 -^]du 

r-.ood( j 8-2)(3- j 9) 7 



1 



oo 

-Kvr/+ , „,\3-/3 „,3-/3 



d(/3-2)(3-/3) 7 



d(/3-2) 
1 



J 

t />oo 



d(/3 



JO 



(* + v) 6 - p - v 6 - p ]dv 
u + v) 2 ~^du dv 

u + v) 2 ~ dvdu. (6.19) 



The result follows on combining (|6.17j) - (|6.19p . □ 

Proof of Theorem [215} Define <p : C R [0,oo) -> C R [0,oo) by 

[<p(x)](t) = X(t) - K I [<p(x)](s)d8. 



Then, (p is a continuous mapping from Ck[0, oo) to Cr[0, oo). 

For any i, T > 0, denote 7£r(t) = ^(T) + Rt(J). Then T&r is a Gaussian process 
with continuous trajectories, and Z(T + •) = <p(TZt)(-)- Therefore, in order to prove that 
Z(T + •) converges in distribution, in Cr[0,oo) to Z^, it suffices to show the convergence 
of TZt to ^oo(O) + cfBh(-), where Zoo(0), and a are as defined before Theorem 12.51 
From Lemmas 16.1116.31 it follows that the finite dimensional distributions of TZt converge 
to those of Z oo (0) + aBH(-). It thus suffices to verify that {TZt(-)}t>o is tight in Cr[0,oo). 
By the Cauchy- Schwartz inequality and (|5.7p . it follows that for any h > 0, t > h and 
T > 0, 

E(\n T (t + h)- n T (t)\\n T (t) - n T (t - h)\) 

= E(\R(T + t + h)- R(T + t)\\R(T + t) - R(T + t-h)\) 

< (E\R(T + t + h)- R(T + t)\ 2 ) 1/2 (E\R(T + t)- R(T + t - h)\ 2 ) 1/2 

~ d(/3-2)(3-/3)(4-/3)- (6 - 20) 
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Note that 4-/3 > 1, since f3 < 3. The desired tightness now follows from standard results 
(cf. Theorems 3.8.6 and 3.8.8 in [2]). This proves the convergence of Z(T + •) to Zoo(-). 
Stationarity of is now immediate. The result follows. □ 



A Auxiliary Results 

Recall that 2 < /3 < 3 and a > /3 — 1. Also recall the notation tf^s) and ^''(r^) 
introduced in the proof of Theorem l2.21 The following lemma gives a key estimate for the 
proof of the theorem. 

Lemma A.l There exists a constant C > depending only on f3 and 6, such that for 
any < h < 1 and h < t < oo, we have the following estimates 

t+h poo 

[^(r,s)] 4 n Q+1 (n0r + l)'^drds < Cn 5 ""- 3 ^ 6 -' 3 , (A.l) 



n^- 1 ) Jo Jo 
and 

>-t+h poo \ 2 



-i / rt+h poo \ * 

4( i_M y [^(r, S )] 2 n Q+1 (n0r+l)-^rd S J < Cn & -^~ 2a h 2 ^ . (A.2) 



Proof: We first prove (jA.ip . We can write the left-hand side of (|A.1|) as 

rt+h poo 



n 4 ( 

1 



-i rt+h poo 

— j t J [r f\(t + h- s)] 4 n a+1 (n9r + l)-Pdrds 

[r A (t + h - s) - r A (i - s)] 4 n a+1 (n9r + l)~^dr ds. (A.3) 



t />oo 



We now bound the two terms in (|A.3|) separately. In the following, we use C > to 
denote a generic constant depending only on j3 and 9; the value of C may change from 
one line to next. Using (n9r + l)"' 3 < (n9r)~@, we have 
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n 4 («-D J t Jo 



t+h poo 

r A (t + h - s)} 4 n a+1 (n9r + l)~ p dr ds 



I ft+h rt+h—s 

-i ft+h roo 
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K 10 ' Jt+h-s 

a-fi+1 rt+h rt+h-s a-(3+l rt+h roo 

< / r^drds+- s -x / (i + fr-*)V*drcfa 



n 4 ^" 1 ) A 7o n 4 («-!) ^ J t+h ^ s 

rt+h 

Cn 5-P-3a / ^ + h _ S f-P ds < £^5-/3-3^6-/3 (A.4) 
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For the second term in (|A.3|) . we have 
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t roo 
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[r A (t + h - s) - r A (t - s)] 4 n a+1 {n9r + l)~^dr ds 



n^' 1 ) Jo Jt-s 

ft roo 



t rt+h—s 

[r-(t- s)] A n a+1 {n6r + l)~ p dr ds 




i rt roo 

H — 77 — tt / / h 4 n a+1 (n9r + l)-Pdrds 
n 4 ^" 1 ) Jo J t+h - s 



1 



t i-h+s 

(r - s) 4 n a+1 (n9r + l)~ p dr ds 




n 4 («-D Jo Js 

-i f't roo 

+— -— T / / frV* +1 (n0r+ir^rds. (A.5) 

For the first term in (|A.5p . we change the order of integration and obtain 

rt rh+s 



— 77- — tt f [ (r-sfn a+l {ner + l)-Pdrds 
n 4 ^" 1 ) Jo Js 



1 



ft 



< A/n _^ I J (r - s) A n a+1 (n9r+-l)' l3 dsdr 
rt+h rr 



+ * / / (r-s) 4 n Q+1 (n0r + l)- /3 ds(ir 
^ Jh Jr—h 

i /-ft i rt+h 

Cn «-^+i /-A Cn a -P +1 f t+h ,,/ 1 \ , 

< — 77 — tt- / r 5_/:, dr H 77 — 77- / h 5 r + — dr 

< Cn 5 -P- 3a h 6 -P, (A.6) 

since {h+^) l ~ P - {t + h+±f~ P < {h+^f' P < h 1 "? for all n G N and 9 6 (0,oo). 
For the second term in (|A.5j) . we have 

1 ft roo /-1 a— rt roo 

1 1 h A n a+l (n9r + l)-^drds < / / h A r~^drds 



= Cn 5 -P- 3a / /i 4 (£ + /t - s) 1 "^ 
J 

< Cn 5 ^- 3a li 6 ^. (A.7) 

Combining <[XT3|) - <[XT7|) . the bound ([ATT]) follows. 

Next we show (|A.2p . The left-hand side of ()A.2p is bounded by 
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+ 



n 4(a-l) 
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4(a-l) 



y y [r/\(t + h-s)] 2 n a+1 (n6r + l)-Pdrdsj 

(rt roo \ 2 

y y [rA(t + /i-s)-rA(i-s)] 2 n a+1 (n0r + l)^drdsJ .(A.8) 
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We bound the first term in (|A.8|) as 

"t+h roo 



4(a-l) 



(pt+h poo \ • 

^ y [rA(t+/i-s)fn a+1 ()ifr + ir^r(lsJ 

4 / pt+h pt+h-s \ 2 
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r .„2(o-/3+l) / /-t+h /-oo n, 2 



t) 4(q-1) x , , 
n v ; V( Jt+h-s 

r-t+h 



< Cn 6 -^' 2a (J (t + h-sf-Pds) <Cn G - 2 ^ 2a h 2 ^\ (A.9) 
For the second term in (\A.8\i , by a change of variables we obtain 

(ft poo \ 2 

y y [rA(t + /i-s)-rA(t-s)] 2 n Q+1 (n0r + l)~ /3 (irdsJ 
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V./n ./,-, 
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4(a-l) 
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4(a-l) 



fe 2 n a+1 (ner+l)~^drds) 

.fft-s / 
pt ph+s \ : 

y y (r-s) 2 n a+1 (n^r + l)- /3 dr(isj 



+ w 4(!-i) (j[ j[ fc 2 n a+1 (ner + l) - ^drds^ . (A.10) 



By changing the order of integration, the first term in (|A,10p is bounded as 

A /pi ph-\-s \ ^ 

\L L ( r - s ) 2na+i ( n9r+i )~ Pdrds ) 

- n 4(«-i) j[ (r~ s) 2 n a+1 (n6r + l)-Pdsdr) 

o / /-t+ft rr \ 2 

+ — —rr / / (r-s) 2 n a+1 (n^r + l)- /3 dsdr 



< Cn 6 " 2 "- 2 ^ 2 ^. (All) 
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The second term in (jA.lOp can be bounded as 

4r) 2(a-/3+l) / r t roc \ 2 / i-t \ 2 

a l, ^-h ^ c " 6 " 2 °"" a " 2{t + " - ^ 

< Cn & - 2a -^h 2 ^. (A.12) 
The bound in (|AT2|) now follows from (|A.8|) - (jA.ll|) . □ 
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